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Abstract. Cosmological equations for homogeneous isotropic models filled by scalar fields and 
• ultrarelativistic matter are investigated in the framework of gauge theories of gravity. Regular 

. inflationary cosmological models of flat, closed and open type with dominating ultrarelativistic 

I matter at a bounce are discussed. It is shown that essential part of inflationary cosmological 

■ models has bouncing character, 
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(N ; 1 Introduction 

^ ■ It is known that inflationary scenario plays the important role in the early Universe theory 
^ ! and allows to resolve a number of problems of standard Friedmann cosmology P3I21- Most 
Q I inflationary models built in the frame of general relativity theory (GR) are singular and 
I limited in the time in the past. In the case of flat and open models such situation is 
^ ■ inevitable 3], and the quantum gravity era cannot be avoided in the past^. At the same 
O ■ time regular inflationary bouncing solutions exist in GR for closed models filled by massive 
O . or nonlinear scalar fields and usual matter [5-12]^. Because limiting (maximum) energy 
density and limiting temperature at a bounce in such models can be essentially less than 
the Planckian ones, classical description of gravitational field in these models is possible 
and quantum gravitational effects are negligible. Note that regular inflationary solutions for 
closed models in GR are unstable at compression stage, and small variations of gravitational 
and scalar fields at compression lead to singularity because of the divergence of the time 
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ci I derivative of scalar field 0. As result the energy density = + V{(j)) and the pressure 
p = |(/)^ — ^(0) (y is a scalar field potential) of scalar field at compression stage near 
singularity are connected in the following way ~ p,/, (Zj. In the case of regular bouncing 
solutions the scalar field potential has to be dominating at a bounce, so in the case of models 
including scalar field and ultrarelativistic matter with the energy density pr the following 
condition at a bounce must be valid 

V - (p^ - Pr > 0. (1) 



^Recently the Pre-Big Bang Scenario in String Cosmology was discussed (See W and Refs. given therein). 
^Eternally inflating regular closed model of an Emergent Universe with asymptotically Einstein static 
state was discussed in Ref. |13| . 
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The relation (1) means that the greatest part of energy density at the end of cosmological 
compression has to be determined by scalar fields^. Because there are not physical reasons 
ensuring the realization of condition (1), regular inflationary solutions in GR cannot be 
considered as a base to build regular inflationary cosmology. 

Different situation takes place in gauge theories of gravitation (GTG) - Poincare GTG, 
metric-affine GTG. Note that GTG are natural generalization of GR by applying the local 
gauge invariance principle to gravitational interaction (see review Jl]). As it was shown in 
a number of our papers [15-20,12], the GTG possess important regularizing properties. By 
satisfying the correspondence principle with GR in the case of usual gravitating systems with 
rather small energy densities and pressures, GTG can lead to essentially different physical 
conclusions in the case of gravitating systems at extreme conditions with extremely high en- 
ergy densities and pressures [13 HI] . In particular, gravitating vacuum with sufficiently high 
energy density p^, > and pressure = —p^ can lead to the vacuum gravitational repulsion 
effect and to a bounce in closed, open and fiat models, that allows to build regular infla- 
tionary cosmological models fHl ^] • Some particular regular inflationary cosmological 
solutions for closed models with scalar fields were discussed in Refs. j211 ITU IT^ . The present 
paper is devoted to study the problem what place take regular inflationary cosmological 
models in GTG. 



2 Generalized cosmological Friedmann equations in GTG 



Homogeneous isotropic models in GTG are described by the following generalized cosmolog- 
ical Friedmann equations (GCFE) 



R^\l-P(p-3p)\ 



3M2 l-(3(p-3p) ' 



(2) 



R + R(lny/\l-(3{p-3p)\ 



An 



(3) 



R 3M2 i-p(p-3p) ' 

where R{t) is the scale factor of Robertson- Walker metrics, k = +1,0,-1 for closed, flat, 
open models respectively, p is energy density, p is pressure, (3 is indefinite parameter with 
inverse dimension of energy density, Mp is Planckian mass. (The system of units with 
h = c = 1 is used). At first the GCFE were deduced in Poincare GTG ^15j, and later it 
was shown that Eqs.(2)-(3) take place also in metric-affine GTG |221I2S!- From Eqs. (2)- (3) 
follows the conservation law in usual form 



P + 3H {p + p) = 0, 



(4) 



•^Note that classical solutions of cosmological Friedmann equations of GR for models including scalar fields 
can be used only for limited time intervals because of the processes of mutual transformations of elementary 
particles and scalar fields in the early Universe. 
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where H = ^ is the Hubble parameter. Besides cosmological equations (2)-(3) gravitational 
equations of GTG lead to the following relation for torsion function S and nonmetricity 
function Q 

S-\Q = -\j^\n\l-(3{p-3p)\. (5) 

In Poincare GTG Q = and Eq. (5) determines the torsion function. In metric-affine GTG 
there are three kinds of models in the Riemann-Cartan space-time {Q = 0), in the 
Weyl space-time {S = 0), in the Weyl-Cartan space-time {S ^ 0, Q ^ 0, the function S 
is proportional to the function Q). The value of determines the scale of extremely 

high energy densities. The GCFE (2)-(3) coincide practically with Friedmann cosmological 
equations of GR if the energy density is small |/9(p — 3p)\ <^ 1. The difference between GR 
and GTG can be essential at extremely high energy densities \f3{p — 3p)\ > 1.^ In the case 
of gravitating vacuum with constant energy density = const > the GCFE (2)-(3) are 
reduced to Friedmann cosmological equations of GR and S = Q = 0, this means that de 
Sitter solutions for metrics with vanishing torsion and nonmetricity are exact solutions of 
GTG j24| ESI and hence inflationary models can be built in the frame of GTG. 

In order to analyze inflationary cosmological models in GTG let us consider systems 
including scalar field minimally coupled with gravitation and usual matter in the form of 
ultrarelativistic matter. This assumption is available by analysis of the hot Universe models 
in the beginning of cosmological expansion. In the case of other form of gravitating matter 
our consideration would be more complicated. If the interaction between scalar field and 
ultrarelativistic matter is negligible, the energy density p and pressure p take the form 



p 



^(0) + gP. 



and the conservation law (4) leads to the scalar field equation 

+ 3E'<\> = -V 



(6) 



(7) 



and the conservation law for matter, which in our case has the following integral prR^ = const. 
By using Eqs. (6)-(7) the GCFE (2)-(3) can be written in the following form 



_k_ 



H 



in 



3M2 



- 3(3V'(j) 



H 



2/9(2^ + 02 



+ 12f3H<i)V' 



Z + H 
2/?(2V^ + 02 
Stt 



1 - Af3{V - 402) 



3/3 



{V"(/)'' - V''')Z + 6/30V'2 



3M2 



1 



(9) 



■^Ultrarelativistic matter with equation of state p = \p exceptional system because Eqs. (2)-(3) are 
identical to Friedmann cosmological equations of GR in this case independently on values of energy density. 
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where Z = 1 - /3UV - 0^), V = V" = ^tt- Relation (5) takes the form 

d(p a(p^ 

S-^Q^2E^_ (10) 
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Unhke GR the cosmological equation (8) leads to essential restrictions on admissible values 
of scalar field and permits to exclude the divergence of derivative (p for any finite value of 0, 
if /3 < 0. Imposing /? < 0, we obtain from Eq. (8) in the case k = 0,+l 

Z>0 or <j)'^ < W + \l3\-\ (11) 

Inequahty (11) is vahd also for open models discussed in Sec. 3. The region E of admissible 
values of scalar field on the plane P with the axis (0, 0) determined by (11) is hmited by 
two bounds L± 

^ ^ ± {AV + (12) 
Prom Eq. (8) the Hubble parameter on the bounds L± is equal to 

H= '^""^ . . (13) 
1 - 2(3{2V + 02) 

According to (13) the right-hand part of Eq. (10) is equal to |if on the bounds L±. This 
means that the torsion (nonmctricity) will be regular, if the Hubble parameter and scalar 
field arc regular. In connection with this our main attention in Sec. 3 will be turned to 
study properties of solutions of GCFE (8)-(9). In accordance with (11) energy density 
and pressure of scalar field satisfy the condition < ^{p^ + \P\~^). Unlike GR the 
equation of state for scalar field p^ ~ p^ is not valid at any stage of models evolution in 
GTG. 



3 Regular inflationary models in GTG 



Let us consider the most important general properties of cosmological solutions of GCFE 
(8)- (9). At first, note by given initial conditions for scalar field (0,0) and values of R and 
Pr there are two different solutions corresponding to two values of the Hubble parameter 
following from Eq. (8): 

3/5\/'0 ± 



H. 



± 



1-2/3(2^ + 02) 



where 



87r 
3MI 



Pr + 7^0 +1^ 



4 



Unlike GR, the values of Hj^ and H_ in GTG are sign-variable and, hence, both solutions 
corresponding to and H_ can describe the expansion as well as the compression in 
dependence on their sign. Below we will call solutions of GCFE corresponding to Hj^ and 
H as i7_|_-solutions and if_-solutions respectively. In points of bounds L± we have D = 0, 
if+ = if_ and the Hubble parameter is determined by (13). The bounds L± are particular 
curves for GCFE. If Z = 0, Eqs. (8)-(9) are satisfied, and because the bounds L± are not 
limited for applying potentials 1^(0) and tend to infinity under \(f)\ — > cxo, corresponding 
solutions of GCFE are singular. 

In order to study the behaviour of cosmological models at the beginning of cosmological 
expansion, let us analyze extreme points for the scale factor R{t): Rq = R{0), Hq = H{0) = 
0. Denoting values of quantities at t = by means of index "0", we obtain from (8)-(9): 



+ gp'v'Ul 



in 



3M2 



Pro + ^0o' + Vo-^(3 (^4Vo - 4^ ' 



(14) 



r stt 

\3MI _ 



Vo-<j>l- Pro - \mVo - 02^2 



+ 3(3 {V\<j)'o-V't)Zo + 6mV", 1 - 2/9(2^0 + 



-1 



^0 5 



(15) 



where Zq = 1 — /3(4Vo — 0o). A bounce point is described by Eq. (14), if the value of Hq is 
positive. By giving concrete form of potential V and choosing values of Rq, (j)Q, 0o and p^o 
at a bounce, we can obtain numerically particular bouncing solutions of GCFE for various 
values of parameter (3. 

The analysis of GCFE shows, that the properties of cosmological solutions depend es- 
sentially on the parameter (3, i.e. on the scale of extremely high energy densities. From 
physical point of view interesting results can be obtained, if the value of \f3\~^ is much less 
than the Planckian energy density fl^, i.e. in the case of large in module values of parameter 
[3 (by imposing Mp = 1). In order to investigate cosmological solutions at the beginning of 
cosmological expansion in this case, let us consider the GCFE by supposing that 

> 1, pr + l<P^ + v<t: m Uv -<p^y . (16) 



puv- 



Note that the second condition (16) does not exclude that ultrarelativistic matter energy 
density can dominate at a bounce pr V + \(t? 12 . We obtain: 

2/7 ( 21/ + 02 ] + 3^' ( 



k_ 
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4V 
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(17) 



( 21/ + 02 ) f 4y _ 02 



( 8V'^ - 341/0^ - 0^ j - 12HV'<j) (2V + 02 



+ 2 ( y"^" 



V 



/2 



4y- 



'2/2 



3M2 



AV - 



(18) 



5 



Eqs. (17)-(18) do not include radiation energy density, which does not have influence on the 
dynamics of inflationary models in the case under consideration (although, as it was noted 
above, the contribution of ultrarelativistic matter to energy density can be essentially greater 
in comparison with scalar field), moreover Eqs. (17)-(18) do not contain the parameter /3. 
According to Eq. (17) the Hubble parameter in considered approximation is equal to 



-3V'd>± 



W- 



2 (2V + 

and extreme points of the scale factor are determined by the following condition 



^0 UVo-<j>l 



2tt 
3MI 



(19) 



Prom Eq. (18) the time derivative of the Hubble parameter at extreme points is 



IT 



3M2 



- -^o) ' + ^ [V^ - V^4) + (aVo - 4) "'1 {2Vo + 4>l) . (20) 



Obviously Eq. (19)-(20) correspond to (14)-(15) in considered approximation. 

Now let us consider flat models, for which Eq. (19) can be written in the following form 



(21) 



Solutions of Eq. (21) arc in physical region S on the plane P determined in considered 
approximation according to (11) by the following inequality 41/ — 0^ > (the neighbourhood 
of origin of coordinates on the plane P is not considered in our approximation). Eq.(21) has 
as solutions two curves Bi and B2 on the plane P, which are situated near bounds L+ and 
L_ respectively. Because of (21) the bounce condition Hq > can be written in the form 



(4Vo-4y + V^-K4>0 or 6(^^y (Fo'4)' + Fo"-n>o>0. (22) 



2tt 

mi 

If the condition (22) is satisfied on the curves Si, 2, the bounce will take place for corre- 
sponding solutions in points of these curves, which can be called "bounce curves". Note 
that bounce curves exist in the case of various scalar field potentials applying in infiationary 
cosmology, in particular, in the case of power potentials V — acf)^'^ {a — const > and n 
is integer positive number). Each of two curves B12 contains two parts corresponding to 
vanishing of or i/_ and denoting by (.81+ , -B2+) and (-Bi-, B2-) respectively. If V' is 
positive (negative) in quadrants 1 and 4 (2 and 3) on the plane P, the bounce will take place 
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in points of bounce curves -Bi+ and B2+ (-Bi- and -B2-) in quadrants 1 and 3 (2 and 4) for 
i7+-solutions (//--solutions) (see Fig. 1). As result practically all phase trajectories for i/+- 
solutions (//--solutions) reaching bounce curves /?i+ and B2+ (/?i- and B2-) correspond 
to bouncing solutions. To obtain regular bouncing solutions we have to take into account 
that bounds L± defined by relation Z = and Eq.(13) correspond to particular solutions of 
GCFE (L-i- -solutions), in points of which //.-solutions reach the bounds and //^.-solutions 
originate from them^. As result the regular transition from //.-solutions to //^.-solutions 
is possible on bounds L±, in particular, in the form of their glueing. In consequence of this 
each solution containing as parts //--solution and //^-solution describes regular inflation- 
ary model, and a bounce takes place by intersection of phase trajectory with corresponding 
bounce curve. However, singular solutions exist also. As it was noted above in the case of 
various potentials applying in inflationary cosmology (in particular, power potentials) the 
curves L± are not limited on the plane P and tend to infinity; the scalar field satisfying the 
following equation on bounds = 2V' diverges in the past and in the future, and hence 
particular L± -solutions are singular. As result any solution containing //_|_-solution (or //_- 
solution) glueed with particular L-|--solutions is singular in the past (or in the future). Note 
the number of such singular solutions is much less than the number of regular solutions, 
namely to one singular solution correspond infinite number of regular solutions obtaining 

^Remind that iJ+ = on the bounds L±. The Hubble parameter is negative in regions between 
curves (L+ and (i- and i?2+), and the value of is positive in regions between curves (L+ and 

Bi-), (i_ and -62-)- In all other admissible regions on the plane P the sign of values iJ+ and iJ_ is normal: 
H+ > 0, iJ_ < 0. 
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by glueing of noted above singular solution with i/^-solution (or if_|_-solution) in different 
points of L±-curve and excluding then divergent part of L± -solution. 

In the case of open and closed models Eq. (19) determines 1-parametric family of bounce 
curves with parameter Rq. Bounce curves of closed models are situated on the plane P 
in region between two bounce curves Bi and B2 of flat models, and in the case of open 
models bounce curves are situated in two regions between the curves: L4. and L_ and 
B2. Because the behaviour of bounce curves for open models is like to that for flat models, 
the situation concerning bouncing inflationary solutions in the case of open models is the 
same as described above situation for fiat models. Unlike fiat and open models, for which 

= only in points of bounds L± and regular inflationary models can be built if 
and i7_-solutions reach boundsL-t, in the case of closed models the regular transition from 
iJ_-solution to i^f+-solution is possible without reaching the bounds L±. It is because by 
certain value of R we have = H_ in the case if Z 7^ 0. Similar regular inflationary 
solution was considered in Ref. 12j. Note that bouncing solutions exist not only in classical 
region, where scalar field potential and kinetic energy density of scalar field go not exceed 
the Planckian energy density, but also in regions, where classical restrictions on scalar fields 
are not fulfilled and according to accepted opinion quantum gravitational effects can be 
essential. 

In general case, when approximation (16) is not valid, bounce curves of cosmological 
models determined by Eq.(14) depend on parameter f3 . By certain value of /3 we have 
1-parametric family of bounce curves with parameter p^o for fiat models, and we have 2- 
parametric families of bounce curves for closed and open models with parameters Rq and 
PrQ. According to Eqs. (14)-(15) the bounce condition Hq > leads to the following relation 



Svr 



Vo + Ipro - \P (4K, - 4) 1-/3 {V^ - KVo) - ^ ^0 > (23) 

It is follows from (23), unlike GR the presence of ultrarelativistic matter does not prevent 
from the bounce realization (compare with (1)) by reaching certain bounce curve. We see the 
GCFE for homogeneous isotropic models including scalar fields and ultrarelativistic matter 
allow to build regular inflationary cosmological models of fiat, open and closed type, although 
singular solutions exist also. The problem of excluding singular solutions in inflationary 
cosmology in the frame of GTG is analyzed in Ref. [26j. and solutions of GTG near the 
origin of coordinates on the plane P are discussed in Ref. j22j . 

As illustration of obtained results we will consider particular bouncing cosmological in- 
flationary solution for fiat model by using scalar field potential in the form V = |m^0^ 
(m = 10~^Mp). The solution was obtained by numerical integration of Eqs. (7), (9) and 
by choosing in accordance with Eq.(14) (or (21)) the following values of scalar field at a 
bounce: 0o = 10^ Mp, 0o = ^3.96757^0 iP = -10^^ M'^). A bouncing solution includes: 
quasi-de-Sitter stage of compression, the stage of transition from compression to expansion, 
quasi-de-Sitter inflationary stage, stage after inflation. The dynamics of the Hubble pa- 
rameter and scalar field is presented for different stages of obtained bouncing solution in 
Figures 2-4 (by choosing Mp = 1). The transition stage from compression to expansion 
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Figure 2: The stage of transition from compression to expansion. 
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Figure 3: Quasi-de-Sitter stage of compression and inflationary stage. 
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Figure 4: The stage after inflation. 

(Fig. 2) is essentially asymmetric with respect to the point t = because of 0o 7^ 0. Similar 
asymmetry is inevitable property of bouncing solutions for flat and open models. In course 
of transition stage the Hubble parameter changes from maximum in module negative value 
at the end of compression stage to maximum positive value at the beginning of expansion 
stage. The scalar field changes linearly at transition stage, the derivative grows at first 
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from positive value 0i ~ 1.6 ■ 10~^ to maximum value ~ 0o and then decreases to negative 
value 4>2 ~ — 1.6 • 10^^. Quasi-de-Sitter inflationary stage and quasi-de-Sitter compression 
stage are presented in Fig. 3. Although the GCFE (8)-(9) and their approximation (17)-(18) 
have different structure from cosmological Friedmann equations of GR, like OR the time 
dependence of functions H[t) and (j){t) at compression and inflationary stages is linear. The 
amplitude and frequency of oscillating scalar field after inflation (Fig. 4) are different than 

< 1 at 
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that of GR, this means that approximation of small energy densities 

the beginning of this stage is not valid; however, the approximation (17)-(18) is not valid 
also because of dependence on parameter [3 of oscillations characteristics, namely, amphtude 
and frequency of scalar field oscillations decrease by increasing of |/?| ^21- The behaviour 
of the Hubble parameter after inflation is also noneinsteinian, at first the Hubble parame- 
ter oscillates near the value H = 0, and later the Hubble parameter becomes positive and 
decreases with the time like in GR. Before quasi-de-Sitter compression stage there are also 
oscillations of the Hubble parameter and scalar field not presented in Figures 2-4. As it 
was noted above^, solutions similar to obtained one can be used only for limited time in- 
tervals. Ultrarelativistic matter, which could dominate at a bounce, has negligible small 
energy density at quasi-de Sitter stages. At the same time the gravitating matter could be 
at compression stage in more realistic bouncing models, and scalar fields could appear only 
at certain stage of cosmological compression. 



4 Conclusion 

As it is shown in the present paper, regular inflationary cosmological solutions for flat, closed 
and open models obtained in the frame of GTG possess some interesting physical properties, 
if the scale of extremely high energy densities is much less than the Plackian one. The 
dynamics of cosmological models at a bounce is determined by scalar fields and Newton's 
gravitational constant; ultrarelativistic matter, which can dominate at a bounce, does not 
have influence on their dynamics. It is shown, that essential part of flat, open and closed 
inflationary cosmological models has bouncing character. 
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